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The effects of pure multiplicative noise on stochastic resonance in an anti-tumor system modulated 
by a seasonal external field are investigated by using theoretical analyses of the generalized potential 
and numerical simulations. For optimally selected values of the multiplicative noise intensity quasi- 
symmetry of two potential minima and stochastic resonance are observed. Theoretical results and 
' numerical simulations are in good quantitative agreement. 
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Chemotherapy remains a traditional option for most advanced cancer. Immunotherapy, however, is a less conven- 
H- 5 tional treatment modality. Usually, chemotherapy and immunotherapy have been regarded as unrelated, so relatively 
^1 , little research has investigated the relationship between these two therapies. Chemotherapy kills tumor cells in a spe- 
■ cial way periodically, but immunotherapy restrains the growth of tumor cells in a more likely linear way [1, 2]. Since 
these different responses of tumor cells to chemotherapy and immunotherapy, when taken together, they imply that 
there is an interesting and significative case for combining chemotherapy and immunotherapy in tumor treatment. 

More than ever, cancer research is now an interdisciplinary effort which requires a basic knowledge of commonly 
used terms, facts, issues, and concepts. In the past decade, many studies have focused on the growth law of tumor 
cells via dynamics approach, specially using noise dynamics [3-9]. Phase transition of tumor growth induced by 
noises is one of the most novel foundations in recent years. Another phenomenon — known as stochastic resonance 
(SR) — shows that adding noise to a system can sometimes improve its ability to transfer information. The basic 
CX three ingredients of stochastic resonance are a threshold, a noise source and a weak input, it is clear that stochastic 
"~ '. resonance is a common case and generic enough to be observable in a large variety of nonlinear dynamical systems 
! ' [10, 11]. Thus, it is reasonable to believe that SR can also occur in a tumor dynamical system. 

The mean field approximate analysis is a conventional theory for SR. It is originally proposed for symmetrical 
' bistable systems with additive noise source [12]. The improvements of the theory of SR have included monostable 
7— I , systems [13, 14], asymmetrical systems [15] and double-noises systems [16-17], but in all these studies the system has 
an additive noise source and an independent external field. To pure multiplicative noise systems, especially to many 
complex dynamical systems, it is still far difficult to solve them exactly, thus numerical methods are comparatively 
convenient options to solve these complex dynamical systems [18, 19]. 

In this letter, chemotherapy and immunotherapy are joined by an anti-tumor model with three elements, which are 
(1) a fluctuation of growth rate, (2) an immune form, and (3) a weak seasonal modulability induced by chemotherapy. 
Based on the analyses on its unique stochastic differential equation and relevant Fokker-Planck equation, we investigate 
a new type of SR phenomenon of anti-tumor model through both theoretical analysis and numerical calculation. We 
1 . call this effect pure multiplicative stochastic resonance (PMSR) to emphasize that pure multiplicative noise induces 
a synchronization, which can be described by the symmetry of the potential wells. A presupposition is given about 
the relationship between SR and tumor response to treatment. 

Lefever and Garay [20] studied the tumor growth under immune surveillance against cancer using enzyme dynamics 
model. The model is, 
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in which X, P, E$ and E are respectively cancer cells, dead cancer cells, immune cells and the compound of cancer 
cells and immune cells, 7, A, ki, &2, £3 are velocity coefficients. This model reveals that normal cells may transform 
into cancer cells, and then the cancer cells reproduce, decline and die out ultimately. This model can be derived to 
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an equivalent single- variable deterministic dynamics differential equation [5], 

^T- = r n x n {\ - ^p~) - tp(x n ) (2) 
dT K n 

where x n is the population of tumor cells; r„ is their linear per capita birth rate and K n is the carrying capacity of 
the environment, respectively. (p(x n ) quantifies the abilities of recognizability and attack which immune cells have to 
tumor cells. Its form is given by <p(x n ) = f3x„/ (t 2 +x„) [5, 21], here (3 is the immune coefficient, e gives a measure of the 
threshold where the immune system is 'switched on'. For convenience, we set x — x n /e, r — r n e, K = K n /e, t = r/e, 
and obtain a non-dimension form of Eq.(2), 

Like most species, seasonal growth is a common feature of tumor cells, especially when they are in a periodic 
chcmothcrapeutic treatment [2]. This means that the growth rate of tumor cells should be a periodic form, for 
example a cosine form. If considering the fluctuation of environment, the growth rate r in Eq.(3) should be rewritten 
as r + A cos(wt) + £(t), where is the Gaussian white noises defined as = 0, = 2M5(t — t'), in 

which M is the noise intensity. The equivalent stochastic differential equation of Eq.(3) will be, 



In the absence of external field, i.e., Aq = 0, if set f(x) = rox(l — x/K) — (3x 2 /(l + x 2 ), and g(x) = x(l — x/K), 
one will obtain the Fokkcr-Planck equation of Eq. (4) , 

dP(x,t) _ d[A(x)P(x,t)} d 2 [B(x)P(x,t)} 

dt dx dx 2 { ' 

in which 

A(x) = f(x) + Mg(x)g(x) 

B(x) = Mg 2 (x) (6) 
The stationary probability distribution of the system is obtained from Eqs.(5) and (6), 

P st{x ) = NeM- 1 ^ 1 ] (7) 
where N is a normalization constant, and the generalized potential 
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The generalized potential U e ff(x) versus the population of tumor cells x is plotted in Fig.l at different noise 
intensities M. Obviously, the potential has two stable states, and its minima are given by A(x) — B (x) = 0, i.e., 
ro(l — x/K) — (3x/(l + x 2 ) — M(l — x/K)(\ — 2x/K) = 0. Here the positions x\, X2 of the potential minima 
are regarded as the inactive state and the active state of tumor cells, respectively. The generalized potential is an 
asymmetric bistable potential wells, and its potentials at x\ and X2 change with the noise intensity M. We define the 
difference of the potentials at X\, x 2 as AU = \U e ff(xi) — U e ff(x2)\ and plot the relationship between AU and the 
noise strength M in Fig. 2, observing a minimum at a nonzero noise level. The figure displays that the symmetry of 
the potential wells is determined by the multiplicative noise intensity. Although they are asymmetry at high and low 
values of noise intensity, the potential wells are quasi-symmetry and AU tends to zero at a suitable noise strength. 
This makes SR possible for the systems with asymmetry potential wells. Due to this unique generalized potential, 
the system modulated by an external field undergoes a special response to multiplicative noise. 

If inputting a seasonal signal A cos(u;t), as shown in Eq.(4), a time series is taken to monitor the anti-tumor 
system response based on a numerical method for stochastic differential equations [19]. At low values of Aq and 
u (i.e., A <C l,to <SC 1), the symmetry of two potential wells decides the synchronization of the probability skips 
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between them mainly. In the presence of additive noise, the relationship between the signal-to-noise ratio (SNR\) 
and the height of the potential barrier (AC/*,) is given by SNRi = AnArk/D, — exp(— 2A.Uk/ D), in which 
Nk is an independent constant, and A, D are external field amplitude and additive noise intensity, respectively [12]. 
This indicates that the height of the potential barrier affects the existence of SR. However, in the absence of additive 
noise, the above formulas of the traditional theory of SR are not suitable for this pure multiplicative noise problem. 
Here the symmetry of the potential has the main effects on the synchronized hopping between the two states, i.e., the 
existence of SR is mainly determined by the symmetry of the potential wells but not the height of the potential barrier. 
Accordingly, for the sake of simplicity, we have proposed an approximate substitute (i?i) for the signal-to-noise ratio 
(SNR), which has the same form as the SNRi mentioned above and is simply written as 

R^NoeM^) (9) 

where TVo is a non-dimension proportionality coefficient. 

Figure 3 shows selected time scries at different multiplicative noise intensities, and their associated power spectral 
intensities are shown in Fig. 4. For a right value of noise intensity, the populations of tumor cells skip back and 
forth between active state and inactive state, indicating the synchronous response of tumor cells to the treatment at 
an optimal fluctuation in their growth rate. The SNR is defined as the ratio of intensity of the peak in the power 
spectral intensity to the height of noisy background at the same frequency. After calculating the SNR, the data are 
combined and presented in Fig. 5. The stochastic resonance induced by multiplicative noise is clear and marked, with 
the maxima of the SNR and R\ approximately at M ~ 0.22. This is a pure multiplicative stochastic resonance. With 
an increase in multiplicative noise intensity, the change trend of SNR matches that of R\ very well. 

In conclusion, we have investigated the stochastic resonance induced by pure multiplicative noise in an anti-tumor 
system. The seasonal factor, reflecting the influence of chemotherapy on tumor cells, is introduced into the conventional 
tumor growth model under immune systems surveillance. On the base of analyses on the asymmetry of the generalized 
potential, we have given a parameter to substitute the SNR, which consists with the numerical results very well. Our 
works offer a method to analyze some complex stochastic differential equations, although they are insufficient to give 
an exact description of real tumor growth. Especially, the synchronous response of tumor cells to chemotherapy is one 
of the novel findings. We expect that these analyses and numerical findings will stimulate theoretical and experimental 
works to verify pure SR in real anti-tumor systems with seasonal treatment. 
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The generalized potential for different intensities of multiplicative noise M = 0.01 (solk 
dot). The remaining parameters are ro = 1.0, j3 = 2.0, and K = 10.0. 
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